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We measure the absolute branching fraction for D — > K~tv using partial reconstruction of 
B° — > D* + X£~ Vfi decays, in which only the charged lepton and the pion from the decay D* + — > D n + 
are used. Based on a data sample of 230 million BB pairs collected at the T(4S) resonance with the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC, we obtain B(D° — > K~-k + ) = 
(4.007 ± 0.037 ± 0.072)%, where the first uncertainty is statistical and the second is systematic. 

PACS numbers: 13.25.Ft, 13.20.He, 13.20.Gd 



The decay D a — » K~ir + [H is a reference mode for 
the measurements of the branching fractions of the D° 
to any other final state. A precise measurement of the 
value of B(D° — > K~ir + ) improves our knowledge of 
most of the decays of the B mesons, and of fundamen- 
tal parameters of the Standard Model. For instance, 
the largest systematic uncertainty on the branching ratio 
B(B° — > D*~l + vi) 1 and the experimental uncertainty on 
the determination of the Cabibbo-Kobayashi-Maskawa 
matrix element V c f, from that semileptonic decay are in- 
duced by the uncertainty on B(D° — > K~ir + ). 

CLEO-c [2| has recently published the most precise re- 
sult on this branching fraction, which is widely used 
We present here a more precise measurement based on a 
different technique. We identify D° — > K~ir + decays in a 
sample of D° mesons from D* + — > D°7r + decays and ob- 
tained with partial reconstruction of B° — > D* + (X)£~Di. 

The data sample used in this analysis consists of an 
integrated luminosity of 210 fb -1 , corresponding to 230 
million BB pairs, collected at the T(AS) resonance (on- 
resonance) and 22fb _1 collected 40MeV below the res- 
onance (off-resonance) by the BABAR detector. The off- 
resonance events are used to subtract the non-BB (con- 
tinuum) background. A simulated sample of BB events 
with integrated luminosity equivalent to approximately 
five times the size of the data sample is used for efficiency 
computation and background studies. 

A detailed description of the BABAR detector is pro- 
vided elsewhere Q. High-momentum particles are re- 
constructed by matching hits in the silicon vertex tracker 
(SVT) with track elements in the drift chamber (DCH). 
Lower momentum tracks, which do not leave signals on 
many wires in the DCH due to the bending induced 
by the 1.5 T solenoid field, are reconstructed solely in 
the SVT. Charged hadron identification is performed by 
combining the measurements of the energy deposition in 
the SVT and in the DCH with the information from 
a Cherenkov detector (DIRC). Electrons are identified 
by the ratio of the energy deposited in the calorime- 
ter (EMC) to the track momentum, the transverse pro- 
file of the shower, the energy loss in the DCH, and the 
Cherenkov angle in the DIRC. Muons are identified in 
the instrumented flux return (IFR), composed of resis- 
tive plate chambers and layers of iron. 

We preselect a sample of hadronic events with at least 
four charged tracks. To reduce continuum background, 
we require that the ratio of the 2 nd to the th order Fox- 
Wolfram [H variables be less than 0.6. We then select a 



sample of partially reconstructed B mesons in the chan- 
nel B° — > D* + (X)£~i/i, by retaining events containing a 
charged lepton (£ = e, (i) and a low momentum pion (soft 
pion, 7T+) which may arise from the decay D* + — > D°irf . 
This sample of events is referred to as the "inclusive sam- 
ple" . The lepton momentum |6J] must be in the range 
1.4 < pg- < 2.3 GeV/c and the soft pion candidate must 
satisfy 60 < p + < 190MeV/c. The lepton and soft pion 
minimum momenta are optimized to minimize uncertain- 
ties due to charm production in B decays and tracking 
errors, respectively. Maximum momentum selections are 
determined by the available phase space. The two tracks 
must be consistent with originating from a common ver- 
tex, constrained to the beam-spot in the plane trans- 
verse to the beam axis. Then we combine pg- , p + and 
the probability from the vertex fit into a likelihood ratio 
variable, optimized to reject BB background. Using con- 
servation of momentum and energy, the invariant mass 
squared of the undetected neutrino is calculated as 

Ml = (£ bcam - E D . - E e f - (p D * + p e ) 2 , (1) 

where i^eam is half the total center-of-mass energy and 
Ei (Ed*) and pi {pd*) are the energy and momentum of 
the lepton (the D* meson). Since the magnitude of the 
B meson momentum, ps, is sufficiently small compared 
to pi and pd*, we set pu — in obtaining Eq. [J As a 
consequence of the limited phase space available in the 
D* + decay, the soft pion is emitted nearly at rest in the 
D* + rest frame. The D* + four-momentum can therefore 
be computed by approximating its direction as that of the 
soft pion, and parameterizing its momentum as a linear 
function of the soft-pion momentum. We select pairs of 
tracks with opposite electric charge for our signal (£ T nf) 
and same-charge pairs (£^ irf) for background studies. 

All events where D* + and l~~ originate from the same 
B-meson, producing a peak near zero in the M 2 dis- 
tribution, are considered as signal candidates. Several 
processes contribute: (a) B° — > D* + £~Di decays (pri- 
mary); (b) B —> D* + (mr)£~~ vt where the D* + (n7r) may or 
may not originate from an excited charm state (_D**) and 
n> 1; (c) B° -> D*+D, D — > i~ X and B° -> D* + t~v t , 
t~ — > £~viv T (cascade); (d) B° — > D* + h~ (fake-lepton) , 
where the hadron (h = 7r, K) is erroneously identified 
as a lepton (in most of the cases, a muon). We also in- 
clude radiative events, where photons with energy above 
1 MeV are emitted by any charged particle using PHO- 
TOS v2.03 [7]. The signal region is Ml > -2 GeV 2 /c 4 
and the sideband is -10 < Ml < -4 GeV 2 /c 4 . 
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The background in the inclusive sample consists of con- 
tinuum and combinatorial BB events, which also include 
events where true D* + and t~ from the two different B 
mesons are combined. We determine the number of signal 
events in our sample with a minimum \ 2 fit to the A'l 2 
distribution in the interval -10 < Ml < 2.5 GeV 2 /c 4 . 
We perform the fit in ten bins of the lepton momentum 
in order to reduce the sensitivity of the result to the de- 
tails of the simulation. In each bin we fix the continuum 
contribution to the off-resonance events, rescaled to ac- 
count for the luminosity ratio between the on- and the off- 
resonance samples, while we vary independently the num- 
ber of signal events from primary, from D**, and from 
combinatorial BB, assuming the shapes predicted by the 
simulation. We fix the contributions from cascade and 
fake-lepton decays, which account for about 3% of the sig- 
nal sample, to the Monte Carlo (MC) prediction. We fit 
eight different sets, divided by lepton kind and run con- 
dition. The reduced % 2 s range between 1.1 and 1.4. Fig- 
ureQJa) shows the result of the fit in the yVf 2 , projection. 
The number of signal events with .M 2 > — 2 GeV 2 /c 4 is 
N incl = (2170.64±3.04(siai)±18.1(sysi))xl0 3 . The sta- 
tistical uncertainty includes the statistical uncertainties 
of the off-resonance and of the simulated events. 
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FIG. 1: The M„ distribution of the inclusive sample, for 
right-charge (a) and wrong-charge (b) samples. The data are 
represented by solid points with uncertainty. The MC fit re- 
sults are overlaid to the data, as explained in the figure. 



We then reconstruct for D° — » K~ir + decays in the in- 
clusive sample. We consider all tracks in the event, aside 
from the t~ and 7T+ , with momenta in the direction trans- 
verse to the beam axis exceeding 0.2 GeV/c. We com- 
bine pairs of tracks with opposite charge, and compute 
the invariant mass Mkk assigning the kaon mass to the 



track with charge opposite the n s charge. The kaon can- 
didate must satisfy a loose kaon identification criterion 
that retains more than 80% of true kaons, while rejecting 
more than 95% of pions. We select events in the mass 
range 1.82 < M Kn < 1.91 GeV/c 2 . We combine each D° 
candidate with the 7r+ and compute the mass difference 
AM = M (K-Tt + K+) - M (K~tt + ). We look for signal in 
the range of 142.4 < AM < 149.9 MeV/c 2 . 

This exclusive sample consists of signal events and of 
the following background sources: continuum, combi- 
natorial BB, uncorrelated peaking D* + and Cabibbo- 
suppressed decays. We subtract the continuum back- 
ground using rescaled off-resonance events selected with 
the same criteria as the on-resonance data. Combinato- 
rial events are due to any combination of three tracks, 
in which at least one does not come from the D* + . We 
determine their number from simulated BB events. We 
normalize the simulated events to the data in the AM 
sideband, 153.5 < AM < 162.5 MeV/c 2 , properly ac- 
counting for the small fraction of signal events (less than 
1%) contained in the sideband. We verify that the back- 
ground shape is properly described in the simulation us- 
ing a sample of Z?* + -depleted events, obtained as follows. 
We use wrong-charge events where the kaon has the same 
charge as the ir s , selected in the A4 2 sideband. More than 
95% of the events so selected in the AM signal region are 
combinatorial background, with a residual peaking com- 
ponent from Cabibbo suppressed decays (K + K~ and tt + 
7r~ , see below). After normalizing the level of the sim- 
ulated events in the sideband, the number of events in 
the signal region is consistent with the data within the 
statistical precision of ±1.3%. 

The background from uncorrelated peaking D* + de- 
cays occurs when the D* + and the l~ originate from 
the two different B mesons. These events exhibit a peak 
in AM but behave as combinatorial background in A4 2 . 
We compute their number in the A4 2 sideband data and 
rescale it to the A4 2 signal region using the A4 2 distri- 
bution of the combinatorial simulated events. 

Cabibbo-suppressed decays D° — > K~K + (tt~7t + ) 
contribute to the peaking background, where one of the 
kaons (pions) is wrongly identified as a pion (kaon). 
Simulation shows that these events peak in AM, while 
they exhibit a broad Mktt distribution. We subtract 
this background source using the simulation prediction. 
It should be noted that the contribution from doubly 
Cabibbo suppressed decays is negligible. Figure [2] shows 
the continuum-subtracted distribution for the data with 
the simulated BB backgrounds overlaid. 

The exclusive selection yields A cxcl = (3.381 ±0.029) x 
10 signal events, where the uncertainty is statistical 
only. The detailed composition of the inclusive and ex- 
clusive data sets is listed in Table HI 

We compute the branching fraction as 

B(D° - K-7T+) = A cxcl /(7V incl Ce ( K^ +) ), (2) 
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FIG. 2: Continuum subtracted AM distribution for data 
(points with error bars) and backgrounds overlaid as ex- 
plained in the figure. 



TABLE I: The inclusive and exclusive samples. 



Source 


Inclusive 


(xlO b ) 


Exclusive 


(xlO 4 ) 


Data 


4.4124 


± 


0.0021 


4.727 


± 


0.022 


Continuum 


0.46 


± 


0.0021 


0.309 


± 


0.017 


Combinatorial BB 


1.7817 


± 


0.0007 


0.819 


± 


0.005 


Peaking 








0.163 


± 


0.008 


Cabibbo-suppressed 








0.055 


± 


0.001 


Signal 


2.1706 


± 


0.0030 


3.381 


± 


0.029 



where E(k-tt+) — (36.96 ± 0.09)% is the D° reconstruc- 
tion efficiency from simulation, and ( ~ 1.033 ± 0.002 is 
the selection bias introduced by the partial reconstruc- 
tion. Only the statistical uncertainties are reported here. 
The bias factor £ accounts for the larger efficiency of the 
inclusive event reconstruction for final states with two 
or fewer tracks from D° decays due to the smaller den- 
sity of hits near the 7r+ track. We study these effects by 
comparing data and simulated distributions of the num- 
ber of charged tracks in each event (ntrk) and of other 
quantities sensitive to the soft pion isolation (angle to 
nearest track and track density within 10° cone around 
the 7T+ direction). We weight simulated events to repro- 
duce the data and recompute the bias. We observe an 
efficiency variation of 0.33% due to n tr k and 0.08% due 
to the other variables. The bias does not depend on some 
other variables (p + , number of 7r+ hits in the SVT). The 
systematic uncertainty due to this selection is ±0.35%. 

The main systematic uncertainty on 7V lncl is due to the 
non-peaking combinatorial BB background. We perform 



the same fit to the I irf background control sample and 
the signal-dominated sample. We take the systematic un- 
certainty in the combinatorial background to be the RMS 
scatter in the ratio, calculated for each bin as shown 
in FigureQJb), °f continuum-subtracted data to the value 
of the combinatorial background determined from the fit, 
resulting in an uncertainty of 0.89%. As first noticed 
in the decays B° -> t~viD+, with D+ -> K*p(aj)ir+, 
constitute a right-charge peaking background, because 
the charged pion is produced almost at rest in the D + rest 
frame. In order to estimate the systematic uncertainty 
due to this peaking combinatorial background, we vary 
its total fraction by ±100% in the BB events in the MC. 
The corresponding systematic uncertainty is ±0.34%. 

We consider systematic uncertainties affecting the sig- 
nal distribution. Final state photon radiation in D° 
decays alters the distribution of M K -^+ and thus affects 
the efficiency computation. We estimate a systematic 
uncertainty of ±0.50% due to the final state photon ra- 
diation by varying by ±30% the fraction of reconstructed 
events in the simulation where at least one photon above 
1 MeV is emitted in the D° — > K~n + decay. 

We also vary, by ±30%, the fractions of cascade and 
fake-lepton decays, which are not determined by the fit. 
Finally, we vary in turn by ±100% the number of events 
from each of the five sources constituting the D** sam- 
ples (two narrow and two broad resonant states, and non- 
resonant D* + ir combinations; these last arc described us- 
ing the model of ref. [9j). We repeat the measurement 
and take the variation as the systematic uncertainty. 

The dominant contribution to the systematic uncer- 
tainty on iV oxcl is due to the charged-track reconstruction 
efficiency. The single charged-track reconstruction effi- 
ciency is determined with 0.50% precision, which corre- 
sponds to ±1.00% overall uncertainty. The efficiency for 
K~ identification is measured with ±0.70% systematic 
uncertainty from a large sample of D* + — > -D°7r+, D a —> 
K~-k + decays, produced in e + e~ — > cc events. To esti- 
mate the systematic uncertainty due to the combinatoric 
background subtraction on _/V oxcl , we first vary the num- 
ber of events from combinatorial background below the 
signal peak by ±1.3%, corresponding to the statistical 
uncertainty obtained from the control sample described 
above. This translates in ±0.3% systematic uncertainty 
on the result. We vary the number of signal events con- 
tained in the sideband by ±30% for background nor- 
malization. This induces a systematic uncertainty of 
±0.16%. We vary the fraction of events from Cabibbo 
suppressed decays by ±10%. The systematic uncertainty 
due to the uncorrelated peaking (from data) is negligible. 

When comparing the simulated Mktx distribution to 
the data in a high purity signal sample (obtained by ask- 
ing, in addition to the other cuts, that the hard pion fails 
K, p and I identification criteria and 0.1435 < AM < 
0.1475 GeV/c 2 ) we observe a slight discrepancy, causing 
±0.56% systematic uncertainty in the reconstruction ef- 
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ficiency. We compute the total relative systematic un- 
certainty of ±1.80% from the quadratic sum of all uncer- 
tainties described above and listed in Table HT1 We cross 



TABLE II: The systematic uncertainties of B(D° -» K~ty + ). 





Source 


5{E)/B (%) 




Selection bias 


±0.35 


jyincl 


Non-peaking combinatorial background 


±0.89 




Peaking combinatorial background 


±0.34 




Soft pion decays in flight 


±0.10 




Fake leptons 


±0.08 




Cascade decays 


±0.08 




Monte Carlo events shape 


in ns 
ztU.Uo 




Continuum background 


±0.05 




D** production 


±0.02 




Photon radiation 


±0.02 


^ycxcl 


Tracking efficiency 


±1.00 




K~ identification 


±0.70 




D° invariant mass 


±0.56 




Photon radiation in D° decay 


±0.50 




Combinatorial background shape 


±0.30 




Combinatorial background normalization 


±0.16 




Soft pion decay 


±0.12 




Cabibbo-suppressed decays 


±0.10 


Total 




±1.80 



check our results using different definitions of the AM 
and M K -^+ signal regions and particle identification. We 
split our data into different sub-samples, depending on 
the run conditions. All the results are consistent. 

In summary, we have measured the absolute branching 
fraction of D° — > K~ir + decay with partial reconstruc- 
tion of B° — » D* + (X)£~I/i, and obtain the result 

B(D° -> K~7T + ) = (4.007 ±0.037 ± 0.072)%, (3) 

where the first uncertainty is statistical and the second 
uncertainty is systematic. This result is comparable in 



precision with the present world average, and it is con- 
sistent with it within two standard deviations. 
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